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Abstract. Supergranules are believed to be an evidence for large-scale subsurface con- 
vection. The vertical component of the supergranular flow field is very hard to measure, 
but it is considered only a few ms _1 in and below the photosphere. Here I present the 
results of the analysis using three-dimensional inversion for time-distance helioseismol- 
ogy that indicate existence of the large- magnitude vertical upflow in the near sub-surface 
layers. Possible issues and consequences of this inference are also discussed. 
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1. Supergranules: Convection cells? 

The photosphere of the Sun is fully covered with supergranules, a convection- 
like cellular pattern. Some properties of supergranules are still being debated 
in literat ure, despite the discovery of this phenomenon more than half cen- 
tury ago ( Hart . 19541 ). Surface properties of supe r granu l es were extensively 



studie d by many authors, e.g.. IWang and Zirinl (I198G ). IStrous and Simon 



(jl998h . lDe~Rosa and Toomrell2004l ). lHirzberger et all (120081 ) and niany_oth 



ers. R ecently, the findings about supergranulation were reviewed by lRieutord and Rincon 
(|2010h . 

The main lesson from dozens of papers published over almost six decades 
is that only the surface properties of supergranules are pretty well known. 
The characteristic cell size is around 30 Mm (the corresponding velocity 
spectrum scale- length is around 38 Mm however), lifetime of supergranules 
is a day or more. Supergranules cannot be observed in white light as the 
temperature within the cell varies by less than 4 K. The plasma flows within 
the supergranules are predominantly horizontal with root-mean-square ve- 
locity ~300 ms _1 , thus supergranular pattern is easily observable in Dopp- 
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lergrams. The vertical (radial, with respect to the spherical Sun) component 
of the supergranular flow is hard to measure, its amplitude is usually within 
measurement error bars. Statistically it has been d etermined that the root - 
mean-square vertical velocity range s from 4 ms" 1 ( Duvall and Birchl . I201CH ) 
to 29 ms" 1 (jHathaway et q/.l . I2OO2I ) . 

The depth structure and origin of the supergranulation is practically un- 



know n. Various local helioseismology methods (for a review see lGizon and Birch 



2005) attempted to study the sub-surface structure of supergranular flow. 
Some studies analysed the flow maps at various depths inferred by helio- 
seismic inversions. They investigated correlation of horizontal divergence of 
the flow at various depths with the surface . T his approach was use d e.g. by 



Duvalll (jl9<ih . IZhao and Kosovichevl (j2003l ) . or lSvanda eiToZl (feoOQ) to name 



a few. Although the numerical values differ for individual studies, a general 
conclusion is that supergranules are conv ection- like cell s exten ding to depths 
of 8-25 Mm with a deep "return flow". IZlebcfk et al\ (1201 ll ) attempted to 
segment supergranules as compact objects in the 4-D space (three for spatial 
coordinates and one for time) from the sequence of helioseismic flow maps. 
The authors tried to use the large sample of 3-D supergranules with their 
histories to capture the evolution of an average supergranule in the time- 
depth domain. These results again support supergranules to be two-layer 
structures loosing the coherence at the depth of some 25 Mm or so. Note 
that all studies cited in this paragraph are based on RLS helioseismic inver- 
sions under ray approximation (geometrical optics) applied to phase-speed 
filtered data and obtained using practically the same code. 

Independent studies, employing ridge- filtered dat a with Bo r n app roxi- 
mation wave ker n els an d OLA inverse methods, e.g. IWoodardl ( 20071 ) and 



Jackiewicz et all (120081 ). did not detect the flow reversal. Moreover, they 
pointed out that any inversion for the flow snapshot deeper than 4-6 Mm 
is dominated by random noise and thus does not reveal any information 
about the de ep supergranular flow . The statement about the noise was also 
confirmed by ISvanda et all (120111 ) using data from numerical simulation of 
solar convection. 

Non-helioseismic methods were also used to draw conc lusions about the 
depth structure of supergranules. E.g., Hathawavl (|2012l ) recently showed 
that supergranules may extend to depths equal to their widths indicating 
that analysis of supergranules with differing sizes may lead to a different 
depth structure. 
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Last but not least: a non-standard helioseismic method was used to in- 
vestigate the vertical component of the supergranular flow ju st under the 

surfac e. A systematic offset in the travel-time shifts measured by lDuvall and Hanasoge 
(|2012l ) for large separations between the measurement points using a spe- 
cial spatio-temporal filtering of the data served as evidence indicating a 
large- amplitude flow under the surface of supergranules. The non-standard 
travel-time measurement procedure was selected to avoid the cross-talk be- 
tween the horizontal and vertical flow in supergranules, which was known 
to be hard to avoid in inverse modelling. Using a simple Gaussian model 
constrained by surface measurements, they predicted a peak in the vertical 
flow of 240 ms- 1 at a depth of 2.3 Mm. 



2. Inverse modelling 



I aimed to verify the results of iDuvall and Hanasogd (120121 ) using a stan- 
dard time-distance inversion pipeline (described and tested in lSvanda et al. 
201 ll ). The pipeline implements Subtr active Optimally Localised Averaging 
(SOLA; Piipers and Thompson , 1992) three-dimensional inversion for time- 
distance (jDuvall et all 119931 ) helioseismology applied to a set of travel-time 
maps. 

The travel times were measured from twelve-hour Do pplergram dat 
acub e s acquired by Heli oseismic and Magnetic Imager (HML lScherrer et al. 
20121 . ISchou et aJl2012l ) in the period of 10 June to 10 July 2011, thus re- 
sulting in 64 consecutive sets of travel-time maps. In each datacube, only 
the disc-centre region (512x512 pixels with pixel size of 1.46 Mm) was 
tracked using the code DRMS_TRACKING (Schunker & Burston, unpub- 
lished) deployed at German Science Center for SDO. The Dopplergram 
datacubes were spatio-temporally filtered with ridge filters, separating /, 
Pi , and P2 mod e s. Th e linearised travel times were measured following 
Gizon and Birchl (12004^ using cen t re-to- annulus and centre-to-quadrant ge- 
ometries (jKosovichev and Duvall , 1997 ) with radii of the annuli 5 to 20 
pixels. 

The SOLA inversion was performed in a Cartesian approximation (x 
and y for the horizontal coordinates and z for the vertical) using sens itiv- 
ity kernels computed in Born approximation (IBirch and Gizonl . 120071 ) for 
geometries consistent with travel-time measurements (ridge-filtered modes 
/ to P2 and centre-to-annulus and centre-to-quadrant averaging with same 
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set of annuli radii). Only the near-surface layers were targeted. Note that 
since only the patch around the disc centre is used for the inverse modelling 
(to prevent effects of foreshortening, which is not accounted for in the sen- 
sitivity kernels), approximation by Cartesian coordinates is justified. The 
vertical velocity then almost coincides with the line-of-sight velocity (the 
angle between the two is less than 24 degrees), except for the oposite sign. 

The SOLA algorithm minimises a cost function consisting of many regu- 
larisation terms. Among those, important terms evaluate the localisation in 
the Sun (an averaging kernel), the level of random noise in the results, and 
cross-talk contributions. By choosing values of the trade-off parameters, one 
can regularise strongly about some of these terms, while relaxing demands 
on the others. 



Due to the random-noise realisations, a normal tomography approach 
cann ot be used to meas ure the vertical component of the supergranular 



flow (ISvanda et all 120111 ), In such case, the signal-to-noise ratio is usually 
in the order of unity, i.e., the resulting maps are very noisy. Even in such 
case, the averaging kernel has many side-lobes which introduce unwanted 
artefacts. The horizontal side-lobes are especially dangerous. Despite being 
"almost negligible" in magnitude, they accumulate a large volume and thus 
represent a significant component of the total volume integral over the whole 
averaging kernel. To avoid such contamination, one has to pick the solution 
with good localisation, no side-lobes, and therefore higher noise level. Such 
estimate of the vertical flow would be completely buried in random noise. 

Luckily, independent flow maps may be considered as independent re- 
alisations of the flow, including the realisation of random noise. When one 
averages over N these independent realisations, one increases the signal-to- 
noise ratio \/jV-times. That is the approach governing this study. It i s fair 
to note that same approach was used also by Duvall and Birch ( 2010l ) and 



Duvall and Hanasogd (120121 ) 



The averaging kernels of the discussed inversion are displayed in Fig. [T] 
It is clear that the flow estimates are averaged mostly over depths of 0- 
2 Mm with some contribution from deeper layers, especially in the case of 
the vertical flow inversion. The centre of gravity of the averaging kernels is 
located at a depth of 1.2 Mm. 
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Figure T. Cuts through all components of averaging kernels from near-subsurface inver- 
sions for horizontal (v x ) and vertical (v z ) flows in average supergranule. Noise levels of 
these inversions are discussed in the main text. Over-plotted contours, which are also 
marked on the colour bar for reference, denote the following: half-maximum of the kernel 
(white) and 5% of the maximum value of the kernel (grey). Components of the averaging 
kernels in the direction of the inversion {K.% and K, z z ) are well localised without horizontal 
side-lobes. It is also evident that the contribution of the cross-talk is negligible as the 
cross-talk averaging kernels (ICy and and K.% and ICy) are practically zeroes. 



3. An average supergranule 

The ensemble averaging approach was implemented to suppress the random- 
noise level and to form the set of travel-time maps representing an 'average 
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supergranule'. The supergranules were stacked about locations of their cen- 
tres. The locations of centres of individual s upergranules were determ ined 
using an algorithm s i milar to that used by Duvall and Birch ( 201dl ) and 
Duvall and Hanasogel (|2012l ). The supergranules were searched for in the 



centre-to-annulus difference travel-time maps measured from /-mode fil- 
tered datacubes for the distance range of 8-12 pixels (0.96-1.44°). After 
smoothing with a Gaussian window with Full-Width-at-Half-Maximum of 
3 Mm to suppress random fluctuations, such maps were sensitive to the 
horizontal divergence of the flow, which is assumed to occur inside super- 
granular cells. 

The pixels with highly negative travel time (which corresponds to large 
divergence signal) were identified in each map. This identification resulted 
a list of potential cell centres. In the next step, the list was cycled and the 
distance to all other potential cell centres was calculated. Should the dis- 
tance between any pair of points be less than 23 Mm, the point with weaker 
signal was removed from the list. In such way duplicities were removed. 

In all 64 datacubes, the procedure identified 5582 supergranule centres. 



3.1. Near-subsurface flows in the average supergranule 

For the inverse modelling, maps with independent realisations of travel times 
were averaged about the location of the supergranular centres. This resulted 
in a single set of travel-time maps, representing travel-time shifts caused 
by average anomalies under the average supergranule. The average travel- 
time maps were inverted for all flow components. The averaging over 5582 
supergranules reduced the noise level by a factor of \/5582 = 76 and led to 
the level of random noise in velocity estimates less than 2 ms _1 . 

The magnitude of the vertical and horizontal flow, as a function of the 
distance R from the centre of the average supergranule, is plotted in Fig. [2] 
The vertical velocity peaks in the cell centre with a value ofll7ms _1 . The 
root-mean-square vertical velocity over the are a of the average superg ranule 



is 21 ms _1 , well within the measurement by lHathawav el al\ ( 20021 ). The 
second peak in the vertical velocity, at a distance of 38 Mm from the cell 
centre, indicates the average location of the neighbouring supergranules in 
my sample and agrees well with the scale-length determined from the ve- 
locity spectrum by other authors. The amplitude of the horizontal velocity 
components (in Fig. [2] only the radial component is plotted) is in agreement 
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with literature. 

The known location of the supergranules allows one to also average the 
measured Dopplergrams and obtain the average line-of-sight velocity, which 
corresp onds to the vertica l veloc ity at the surface. My results are similar to 
those of Duvall and Birch i.e. 12±1 ms 1 in the cell centre, and are 

also plotted in Fig. [2] 
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Figure 2: Azimuthal average of the radial (with respect to the cell centre, diamonds) 
and vertical (v z , squares) near-subsurface velocity in the average supergranule. The line- 
of-sight velocity («j os with negative sign; full circles) obtained from the stacked Dopp- 
lergrams of supergranules is over-plo tted for reference. It is practically identical to that 
presented in lDuvall and Birchl l|201(t l. Error bars for all plots are 2 ms -1 or less. 



3.2. IS THIS MODEL OF AN AVERAGE SUPERGRANULE REASONABLE? 

The large-magnitude vertical flow resulting from the inversi on is somewhat 

surpri sing, as this has not been seen before except for study of lDuvall and Hanasoge 
(|2012[ ). Thus, one can ask the question, whether such solution to the prob- 
lem is physically feasible. The check whether the mass conservation holds 
is one of the simplest. By splitting the continuity equation 

V-pv = 0, (1) 
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where p stands for density and v for the full velocity vector, into the hori- 
zontal and vertical parts, one obtains 

P ( Z )V h .v h + ^ = 0. (2) 

To check whether the mass conservation holds, one can use the ratio of 
the two terms, which should be close to unity, and both terms should have 
opposite signs. Here density p is assumed to vary only with depth z, hor- 
izontal variations are neglected. Operator is the horizontal part of the 
divergence operator applied to the horizontal velocity vector = (v x ,v y ) 
and can be evaluated accurately in the wave-vector space. 

d h (x, y) = V h ■ v h (x, y) = J 7 ' 1 [k ■ v h (k x , k y )} , (3) 

where J 7-1 represents the inverse Fourier transform, k = (k x ,k y ) is a hori- 
zontal wave-vector and Vh a Fourier image of the horizontal velocity vector. 

The vertical derivative term can be roughly estimated from the values 
measured at the target depth and at the surface. The extent to which the 
mass conservation holds is evaluated from the equation 

\p(zo)vz(zo) ~ p{zt)vz{zt)]/{z t - go) , . 

p(z t )d h 

where Zt = 1.2 Mm is a target depth and zq = Mm at the surface. In the 
centre of average supergranule, i.e. at (x,y) = (0,0), for t he vertical compo 



nent of the flow one uses value of v z (Q, 0; Zq) = 10 ms 1 ( Duvall and Birch 



2010 ) , which results from direct line-of-sight velocity measurements and thus 



is representative for vertical velocity at the surface. Lets take v z (0, 0;zt) = 
117 ms -1 at the t arget depth and rf/,,(0, 0) = 83.9 m/ s/Mm. Densities are 
given by Model S (jChristensen-Dalsgaard et all Il996l ) . The ratio \x should 
be close to unity if the mass conservation roughly holds. Indeed, in the case 
of here presented model, fj, = 1.16. The vertical derivative of the density mo- 
mentum is only 16% larger than the "compensating" horizontal divergence of 
the density momentum. The best agreement is achieved when one assumes 
the target depth of 1.4 Mm, which, given the rather extended averaging 
kernel, is still a representative value. 

The situation for the off-centre points is more difficult for two reasons. 

1. The averaging kernels for horizontal and vertical components are not 
identical, thus the inversion for independent components cannot be 
compared in a straightforward way. 
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Figure 3: An azimuthal average of the mass-conservation quantity fi defined by Eq. @ 
with 1-cr error-bars as a function of the distance from the average cell centre R. In the 
ideal case, /j, = 1 everywhere (dashed line). 



2. The representative of the surface vertical flow is taken directly from 
the stacked Dopplergrams. The difficulty with the latter can directly 

be seen in Fig. [2] v\ os turns negative much closer to the centre of 

the cell than v z {zt). 

The plot in Fig. [3] shows that |/i| < 10 within the full cell with a change of 
sign, which is correlated to the change of sign of the horizontal divergence 
of the density momentum. Given the roughness of the estimate and taking 
measurement error-bars into account, the mass conservation holds well. 



4. Concluding remarks 

The presented results can still be affected by various biases discussed below. 

• Detection of supergranules The searching algorithm used in this 
study prefers the supergranules with highest divergence signals. It de- 
tected 5582 cells. Assuming that supergranules cover the whole Sun 
in a hexagonal pattern with a flat-to-flat distance of 30 Mm, around 
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50 000 cells are expected. Thus, to obtain a better representative of 
average supergranule a different segmentation algorithm (such as by 



Hirzberger et all 120081 ) must be used, which detects not only the cen- 
tres, but the boundaries as well. It also is not clear whether the large 
divergence signal belongs to supergranulation or to some smaller scale 
(and thus larger amplitude) convection mode. The algorithm based 
on pattern detection would overcome this issue too. 

Stacking of supergranules The detected supergranules are stacked 
to the point of the maximal divergence, which does not necessary 
need to be in their geometrical centre. To overcome this issue, the 
detection code was modified to fit paraboloids in the vicinity of points 
with largest negative difference /-mode travel times and to choose 
vertices of these paraboloids as locations of the cell centres. The results 
remained largely unchanged. 

Travel-time measurement problem Due to the complexity of the 
procedure, it is possible to determine travel times incorrectly, e.g., 
when one makes an error in the implementation of one of the steps nec- 



essary in the travel-time measurement. However. iDuvall and Hanasoge 



(2012F) used a completely different travel-time code and a very differ- 



ent spatio-temporal filtering (phase-speed filters instead of ridge filters 
used here). Despite these significant differences in the implementation, 
they arrived basically at the same results. 

Problem in the inversion code The in yerse problem fo r time - 



distance helioseismology is well described (jJackiewicz et all 120121 ) . 
but an error can be made in the implementation. It has to be noted 
that the same inversi on code was properl y and successfully validated 



using synthetic data (jSvanda et all l201ll ). 



Problem with sensitivity kernels It is possible that the sensitiv- 
ity kernels used in this study resulting from the forward modelling are 
not good enough, which would in turn af f ect th e results of the inverse 



modelling. Again, IDuvall and Hanasogd (120121 ) did not perform the 



inverse modelling and thus did not use any sensitivity kernels in inter- 
preting their results. And again, despite the crucial deviations in the 
approach, the results of their study are qualitatively consistent with 
mine. 
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There is a lot to be done to confirm the mo del of supergranule with large- 
amplitude flows. iDuvall and Hanasogd (120121 ) predicted also large-amplitude 
horizontal flows (some 700 ms" 1 at the depth of 1.6 Mm). Only a very 
detailed tomography with a fine depth resolution may answer related ques- 
tions. 
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